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Abstract 

We investigate the low-energy behavior of the four-point Green's function 
T^ v describing virtual Compton scattering off the nucleon. Using Lorentz 
invariance, gauge invariance, and crossing symmetry, we derive the leading 
terms of an expansion of the operator in the four-momenta q and q' of the ini- 
tial and final photon, respectively. The model-independent result is expressed 
in terms of the electromagnetic form factors of the free nucleon, i.e., on-shell 
information which one obtains from electron- nucleon scattering experiments. 
Model-dependent terms appear in the operator at O(q a q'o), whereas the orders 
0(q a qp) and 0(q' a q'g) are contained in the low-energy theorem for , i.e., no 
new parameters appear. We discuss the leading terms of the matrix element 
and comment on the use of on-shell equivalent electromagnetic vertices in the 
calculation of "Born terms" for virtual Compton scattering. 
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I. INTRODUCTION 



Low-energy theorems (LET) play an important role in studies of properties of parti- 
cles. Based on a few general principles, they determine the leading terms of the low-energy 
amplitude for a given reaction in terms of global, model-independent properties of the par- 
ticles. Clearly, this provides a constraint for models or theories of hadron structure: unless 
they violate these general principles they must reproduce the predictions of the low-energy 
theorem. On the other hand, the low-energy theorems also provide useful constraints for 
experiments. Experimental studies designed to investigate particle properties beyond the 
global quantities and to distinguish between different models must be carried out with suf- 
ficient accuracy at low energies to be sensitive to the higher-order terms not predicted by 
the theorems. Another option is, of course, to go to an energy regime where the low-energy 
theorems do not apply anymore and model-dependent terms in the theoretical predictions 
are important. 

The best-known low-energy theorem for electromagnetic interactions is the theorem for 
"Compton scattering" (CS) of real photons off a nucleon [Q-^]. Based on the requirement of 
gauge invariance, Lorentz invariance, and crossing symmetry, it specifies the terms in the low- 
energy scattering amplitude up to and including terms linear in the photon momentum. The 
coefficients of this expansion are expressed in terms of global properties of the nucleon: its 
mass, charge and magnetic moment. In experiments, one can make the photon momentum 
- the kinematical variable in which one expands - small to ensure the convergence of the 
expansion and to allow for a direct comparison with the data. By increasing the energy of 
the photon one will become sensitive to terms that depend on details of the structure of the 
nucleon beyond the global properties. Terms of second order in the frequency, which are not 
determined by this theorem, can be parameterized in terms of two new structure constants, 
the electric and magnetic polarizabilities of the nucleon (see, for example, Q). 

As in all studies with electromagnetic probes, the possibilities to investigate the structure 
of the target are much greater if virtual photons are used. A virtual photon allows one to 
vary the three-momentum and energy transfer to the target independently. Therefore it 
has recently been proposed to also use "virtual Compton scattering" (VCS) as a means to 
study the structure of the nucleon [§-0]. The proposed reaction is p(e, e'p)j, i.e., in addition 
to the scattered electron also the recoiling proton is detected to completely determine the 
kinematics of the final state consisting of a real photon and a proton. It is the purpose of 
this note to extend the standard low-energy theorem for Compton scattering of real photons 
to the general case where one or both photons are virtual. The latter would be the case, e.g., 
in the reaction e~ + p — > e~ + (e~ + e + ) + p. We will refer to both possibilities as "VCS" . 

There are several different approaches to derive the LET for Compton scattering of real 
photons. One was first used by Low |]. It made use of the fact that in terms of "unitarity 
diagrams" the scattering amplitude is dominated by the single-nucleon intermediate state. 
In such unitarity diagrams - not to be confused with Feynman diagrams - all intermediate 
states are on their mass shell M. Lorentz invariance and gauge invariance then allow a 
prediction for the amplitude to first order in the frequency. Another approach , first used 
by Gell-Mann and Goldberger ||, relies on a completely covariant description in terms of the 
basic building blocks of electromagnetic vertices and nucleon propagators. They split the 
amplitude into two classes, A and B: class A consists of one-particle-reducible contributions 
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that can be built up from dressed photon-nucleon vertices and dressed nucleon propagators. 
Class B contains all one-particle-irreducible two-photon diagrams, where the second photon 
couples into the dressed vertex of the first one. Application of two Ward-Takahashi identities, 
one relating the photon-nucleon vertex operator to the nucleon propagator, the other relating 
the irreducible two-photon vertex to the dressed one-photon vertex, then lead to the same 
result for the leading powers of the low-energy amplitude. One can use yet another technique 
introduced by Low M to describe bremsstrahlung processes. This technique relies on the 
observation that poles in the photon momentum can only be due to photon emission from 
external nucleon lines of the scattering amplitude. Low's method has, in a modified form, 
also been widely used in the framework of PCAC |L0] . 

So far, there have been only a few investigations of the general VCS matrix element, 
since most calculations were restricted to the Compton scattering of real photons. In JTTJ 
electron-proton bremsstrahlung was calculated in first-order Born approximation. The pho- 
ton scattering amplitude, for one photon virtual and the other one real, was analyzed in 



terms of 12 invariant functions of three scalar kinematical variables. In [|T^] it was shown 
that the general VCS matrix element - virtual photon to virtual photon - requires 18 invari- 
ant amplitudes depending on four scalar variables. In [13|] the reaction 7p — > p + e + e~ was 
investigated. Sizeable effects on the dilepton spectrum from the timelike electromagnetic 
form factors of the proton were found. Very recently the low-energy behavior of the VCS 



matrix element was investigated ||14|| . Using Low's approach |9j the leading terms in the 
outgoing photon momentum were derived. It was shown that the virtual Compton scatter- 
ing amplitude at low energies involves 10 "generalized polarizabilities" which depend on the 
absolute value of the three-momentum of the virtual photon. These new polarizabilities were 
estimated in a nonrelativistic quark model. In [15] the VCS amplitude was calculated in the 
framework of a phenomenological Lagrangian including baryon resonances in the s and u 
channel as well as tt° and a exchange in the t channel. A prediction for the \q\ dependence 
of the electromagnetic polarizabilities a and f3 was made. 

The purpose of this work is to identify, in an analogous form to the real CS case, those 
terms of VCS which are determined on the basis of only gauge invariance, Lorentz invari- 
ance, crossing symmetry, and the discrete symmetries. In the following, we will refer to 
such terms not fixed by this LET for simplicity as "model dependent". By introducing 
additional constraints, such as chiral symmetry, also statements about these terms become 
possible. This is, however, beyond the scope of the present work. In our study of low-energy 
virtual Compton scattering, below the onset of pion production, we will mainly work on 
the operator level. This allows us to work without specifying a particular Lorentz frame or 
a gauge. We combine the method of Gell-Mann and Goldberger [EJ with an effective La- 
grangian approach. Class A is obtained in the framework of a general effective Lagrangian 
describing the interaction of a single nucleon with the electromagnetic field |@. In the spe- 
cific representation we choose, this turns out to be a simple covariant and gauge invariant 
"modified Born term" expression, involving on-shell Dirac and Pauli nucleon form factors Fi 
and F 2 . The Ward-Takahashi identities then allow us to determine the leading-order term 
of the unknown class B contribution in an expansion in both the initial and final photon 
momenta. Furthermore, with the help of crossing symmetry a definite prediction can be 
made concerning the order at which one expects model-dependent terms. 

Our work is organized as follows. We start out in Sec. II by outlining the general structure 
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of the VCS Green's function in the framework of Gell-Mann and Goldberger. We state the 
ingredients for the derivation of the LET, namely, crossing symmetry and gauge invariance. 
Section III derives the LET for virtual photon Compton scattering and we discuss the leading 
terms of the matrix element for the reaction e~ +p — > e~ + p + r y in the center-of-mass frame. 
As the notion of "Born terms" is important for the LET, we comment on this aspect in Sec. 
IV and point out ambiguities that arise in their definition, both for real and virtual photons. 
Our results are summarized and put into perspective in Sec. V. 



II. STRUCTURE OF THE VIRTUAL COMPTON SCATTERING TENSOR AND 

GAUGE INVARIANCE 

In this section we will define the Green's functions and the kinematical variables relevant 
for the discussion of VCS off the proton. We will consider the constraints imposed by the 
fundamental requirements of gauge invariance, Lorentz invariance and crossing symmetry. 
We do this in the framework of a manifestly covariant description, incorporating gauge 



invariance in its strong version, namely, in the form of the Ward-Takahashi identities |16l , |T7 
The approach is similar to that of [^] using, however, a somewhat more modern formulation. 
The electromagnetic three-point and four-point Green's functions are defined as 

G» p {x,y,z) = <0\T(* a (x)*p(y)J»{z))\0>, (2.1) 
G%{w,x,y,z) = <0\T(^ a (w)^ p (x)J^(y)r(z))\0>, (2.2) 

where J M is the electromagnetic current operator in units of the elementary charge, e > 0, 
e 2 /47r = 1/137, and where ^ denotes a renormalized interpolating field of the proton; 
T denotes the covariant time-ordered product [IE]. Electromagnetic current conservation, 



<9 M J M = 0, and the equal-time commutation relation of the charge density operator with the 
proton field, 

[J°(x), *(y)]6(x° - y°) = -5\x - y)*{y), (2.3) 



are the basic ingredients for deriving Ward-Takahashi identities [nj,[T7]. 

Using translation invariance, the momentum-space Green's functions corresponding to 
Eqs. ( |2.1| ) and ( |2.2| ) are defined through a Fourier transformation, 

(2n) 4 5 4 (p f -p l -q)G^(p f ,p l ) = J d'xd'ydhe^-Prv-^G^y^^ (2 . 4) 

(2vr) 4 5 4 ( P/ + q'- Pl - q)G^(P, q, q') = J d 4 wd 4 xd 4 yd 4 ze l ^- w -^ x -^'^G^(w, x, y, z), 

(2.5) 

where pi and pf refer to the four- momenta of the initial and final proton lines, respectively, 
P = Pi + Pf, and where q and —q' denote the momentum transferred by the currents J M and 
J u , respectively. We note that G^p depends on two independent four-momenta, e.g., Pi and 
Pf. In particular, it is not assumed that these momenta obey the mass-shell condition p\ = 
p 2 f = M 2 . Similarly, G^p depends on three four-momenta which are completely independent 
as long as one considers the general off-mass-shell case. This will prove to be an important 
ingredient below when analyzing the general structure of the VCS tensor. 
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Finally, the truncated three-point and four-point Green's functions relevant for our dis- 
cussion of VCS are obtained by multiplying the external proton lines by the inverse of the 
corresponding full (renormalized) propagators, 

r»(p f , Pi ) = [iS{ Vf )}- l G^ Vf , Vi )[iS{ Vl )}-\ (2.6) 
r^(P,q,q') = [iS( Pf )]- l G^{P )qi q')[iS( Pl )]-\ (2.7) 

where, for convenience, from now on we omit spinor indices. Using the definitions above, it 
is straightforward to obtain the Ward-Takahashi identities 

q^(p f ,p i ) = S- 1 (p f )-S- 1 (p i ), (2.8) 
q^(P, q, q>) = i (S-\p f )S(p f - q)V(p f - q, Pi ) - T u (p f , Pi + q)S(p t + q^faj) . (2.9) 

Following Gell-Mann and Goldberger [[J, we divide the contributions to T^ u into two 
classes, A and B, T^ u = +Tg > , where class A consists of the s- and u-channel pole terms; 
class B contains all the other contributions. We emphasize that this procedure does not 
restrict the generality of the approach. The separation into the two classes is such that all 
terms which are irregular for g M — > (or g' M — » 0) are contained in class A, whereas class 
B is regular in this limit. Strictly speaking, one also assumes that there are no massless 
particles in the theory which could make a low-energy expansion in terms of kinematical 
variables impossible furthermore, the contribution due to t-channel exchanges, such as 
a 7T°, has not been considered. 

The contribution from class A, expressed in terms of the full renormalized propagator 
and the irreducible electromagnetic vertices, reads 

IT = r(p />P/ + fliSfa + q)T^ Pi + q, Pi ) + r»(p fiPf - q)iS(pi - Wfa - q', Pi ). 

(2.10) 

Note that is symmetric under crossing, q <-» —q' and \i <-> v, i.e., 

T^(P,q,q') = r/(P,-q',-q). (2.11) 

Since also the total is crossing symmetric, this must also be true for the contribution 
of class B separately @. Using the Ward-Takahashi identity, Eq. (|2.8|) , one obtains the 
following constraint for class A as imposed by gauge invariance: 

q^(P, q, q') = i (T»(p f ,p f + q') - T»{ Pi - q', Pi ) 

+S~ 1 (p f )S(p t - q')T v ( Pl - q', Pi ) - V v { Pf:P} + q')S( Pl + q)^ 1 ^)) 

= r A (P,q,q'). (2.12) 

Similarly, contraction of with q' y results in 

q' u r%{P, q, q') = -i (T»(p f ,p f - q) - T»( Pl + q, Pi ) 

+S- 1 (p f )S(p i + q)T»( Pl + q, Pi ) - r»(p f ,p f - q)S(pi - gO^fo)) 
= -fi(P,-<f,-q), (2-13) 
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which is, of course, the same constraint which one obtains from Eq. ( [2.12| ) using the crossing- 
symmetry property of 

qlr%(P,q,q') = -(-q' v )r/(P,-q',-q) = -f»(p,-q>,-q). (2.14) 



Combining Eqs. (|2.9|) and Q2.12|) generates the following constraint for the contribution of 
class B: 

q„T% = q^ - IT) = i (T»( Pl - q', Pi ) - T^p f ,p f + q')) , (2.15) 

relating it to the one-photon vertex || . Once again, the second gauge- invariance condition, 
obtained by contracting with q' u , is automatically satisfied due to crossing symmetry. 

The 4x4 matrix of class B is a function of the three independent four-momenta 
P, q, and q'. It is important to realize that the twelve components of these momenta are 
independent variables only for the complete off-shell case, i.e., if one allows for arbitrary 
values of p 2 and p 2 . This will be important when making use of the constraints imposed by 
gauge invariance. Using Lorentz invariance, gauge invariance, crossing symmetry, parity and 
time-reversal invariance, it was shown in |T2[ that the general T^ u for VCS off a free nucleon 
with both photons virtual consists of 18 independent operator structures. The functions 
associated with each operator depend on four Lorentz scalars, e.g., q 2 , q' 2 , v = P ■ q = P ■ q', 
and t = (pi — pf) 2 . When allowing the external nucleon lines to be off their mass shell, one 



will have an even more complicated structure [|19 



However, in our derivation of the low-energy behavior of the electromagnetic four-point 



Green's function we will not require the full structure as discussed in JT^]. At low energies 
we expand in terms of the four- momenta q^ and q ,fl , 

T% = a^{P) + b^ p {P)q p + c^{P)q' a + ■■■, (2.16) 

where the coefficients are 4x4 matrices and can be expressed in terms of the 16 independent 



Dirac matrices 1, 75, 7^, 7^75, a^". An expansion of the type of Eq. ( 2.16|) is expected to work 
below the lowest relevant particle-production threshold, in this case the pion-production 
threshold; we refer the reader to, e.g., ]2D| where a similar discussion for the case of pion 
photoproduction can be found. 

So far we have considered general features of the operators entering into the description 
of VCS. It is clearly one of the advantages of using a covariant description of the type of Eq. 



(^■16j) that it neither uses a particular Lorentz system nor a specific gauge. When referring 



to powers of q or q', we mean the ones coming from the Dirac structures and their associated 
functions. This is different when one works on the level of nucleon matrix elements or the 
invariant amplitude, where also kinematical variables from the spinors or normalization 
factors enter into the power counting. 

We conclude this section by noting that the above framework can easily be applied to 
VCS off a spin-0 particle, such as the pion. In that case P*", of course, has no complicated 
spinor structure. The building blocks for class A are simply the corresponding irreducible, 
renormalized electromagnetic vertex for a spinless particle and the full, renormalized prop- 
agator A(p) |2J. 
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III. LOW-ENERGY BEHAVIOR OF VCS 



In a two-step reaction on a single nucleon, such as 7*iV — > 7*iV, the intermediate nucleon 
lines in the s- and u-channel pole diagrams of class A are off mass shell, while the external 
nucleons are on shell. The early, manifestly covariant derivations of the low-energy theorem 
for Compton scattering took into account that the associated half-off-shell electromag- 
netic vertex of the nucleon has a different and more complicated structure than the free 
vertex. However, it was shown that the model- and representation-dependent properties of 
an off-shell nucleon do not enter in the leading terms of the full Compton scattering ampli- 
tude when the irreducible two-photon amplitude of class B is included consistently. This 
was explicitly shown by expanding Eq. (|2.10|) to first order in q^ and by constructing the 
leading-order term of Y^ with the help of Eq. ( [2. 15 ) and crossing symmetry or, equivalently, 
the second gauge-invariance constraint. The final result for the amplitude was given in the 
laboratory frame and in Coulomb gauge. 

In order to obtain the LET for VCS, we will proceed on the operator level and combine 
the method of 0|| with ideas of an effective Lagrangian approach to Compton scattering 
H. Let us first recall that the electromagnetic three-point and four-point Green's functions 
of Eqs. (2.1) and (|2.2 ) depend on the choice of the interpolating field \& of the proton, i.e., 
are "representation dependent". Therefore, the truncated Green's functions in momentum 
space, Eqs. ( |2.6| ) and (|2.7|), are in general not directly related to observables, except forp^ = 
p 2 = M 2 . Consequently, the separation into class A and class B is necessarily representation 
dependent, since the total Green's function, Eq. (p^6|), as well as the individual building 
blocks of class A are representation dependent; this has of course no effect on the final on- 
shell result, which is representation independent. This was explicitly shown in |2^] for the 
case of real Compton scattering off the pion. On the other hand, for any given appropriate 
interpolating field satisfying the equal-time commutation relation of Eq. Q2.3| ) the Ward- 
Takahashi identities, Eqs. ( |2.8| ) and (|2.9| ), hold, providing important consistency relations 
between the different Green's functions. In the following we will make use of these relations 
for arbitrary P, q, and q', which, in particular, includes arbitrary pf and p 2 . Only at the 
end, the observable on-shell case will be considered. 



A. Derivation of the Low-Energy Theorem 

We will derive this theorem by using a convenient representation, the "canonical form" , 
of the most general and gauge invariant effective Lagrangian. In the framework of effective 
Lagrangians, a canonical form is defined as a representation with the minimal number of 
independent structures (see, e.g., p3| ). For the class A terms, we need the electromagnetic 
vertex and the propagator of the nucleon. Below the pion-pro duct ion threshold, the La- 
grangian for a single proton, interacting with an electromagnetic field, can be brought into 
the canonical form 0] 

oo 

C lNN = #(z#- M)tf - eY J {{-^) n ' 1 d v F tlv )F l j! 1 ^ 

n=l 

„ / oo \ _ 

^ + E((- D ) n ^)^2n W*, (3.1) 



AM 



n=l 
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where = (c^ + ieA^)^/, F^ = d ii A v — dyA^. The electromagnetic structure of the 

proton is accounted for through the Dirac and Pauli form factors, Fx and F 2 , respectively, 
which are expanded according to 

oo oo 

F 1 (q 2 ) = l + J2(q 2 ) n F ln , F 2 (q 2 ) = K+J2(q 2 ) n F2n, k = 1.79. (3.2) 

n=l n=l 

To this representation of the Lagrangian belongs, of course, a particular canonical La- 
grangian of order e 2 that generates the class B terms. However, as will be seen below, we 
do not need to know it in detail for this derivation. Clearly, Eq. (|3.1j ) is invariant under the 
gauge transformation ^ — * exp(— iea(x))ty, and A^ — > A^ + d^a. In order to arrive at Eq. 
( pTTj) , use has implicitly been made of the method of field transformations (see, e.g., p3|-^|). 
It should be stressed that all ingredients needed for Eq. ( |3.1| ) are on-shell quantities that can 
be determined model independently from electron-proton scattering. Any explicit off-shell 
dependence of the irreducible three-point Green's function has been transformed away and 
will thus not show up in the class A contribution. Such transformations, however, generate 
concomittant irreducible class B terms for the amplitude that must be treated consistently 
(see [§J for details). 

Using standard Feynman rules, the irreducible vertex associated with the effective La- 
grangian of Eq. ( |3.1| ) is found to be 

q 2 

(3.3) 



r&/(p/,ft) = rXr/(p/ - Pi) = YF.iq 2 ) + pVrf + i w F2(g2) ' q = Pf ~ p 



Since it is an important ingredient in our derivation of the LET, we emphasize that the 
vertex of Eq. ( |3.3| ) satisfies the Ward-Takahashi identity of Eq. Q2.8| ); the corresponding 
propagator in the representation that yields Eq. (|3.1| ) is the Feynman propagator of a point 
proton. As a consequence of gauge invariance, Eq. ( |3.1| ) automatically generates a term in 
the vertex proportional to q^. 

We note that the effective Lagrangian approach provides a natural explanation for the 
vertex of Eq. ( |3.3j ) which has previously been used by several authors as a simple means 
to restore gauge invariance in the form of the Ward-Takahashi identity (see, for example, 
|27|). However, it is not an independent building block for any amplitude, but must be used 



together with the corresponding irreducible class B terms for the reaction in question! 

In principle, we could now proceed to construct the most general effective Lagrangian 
relevant to generate the corresponding class B terms for VCS. This would allow us, together 
with a calculation of the pole terms involving the vertex of Eq. ( ^.3[ ), to determine the model- 
dependent terms of T^ u at low energies. However, at this point it is more straightforward to 
apply the results of the last section. We obtain for in the framework of Eqs. ( [3.1|) and 



IX// = Kff{-q')iS F { Pi + q)K ff (q) + r^ f {q)iS F (pi - q)K ff (-q'). (3.4) 
Making use of Eq. (|2.15|) to obtain the gauge-invariance constraint for class B in our repre- 
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sentation, we see that it has the simple formF] 

q^ = 0. (3.5) 

This is due to the fact that the vertex of Eq. ( |3.3j ) depends on the momentum transfer, only. 
Note that Eq. ( |3.5|) is still an operator equation, valid for arbitrary P, q, and q' . 

Class A is defined to contain for any representation the pole terms of T^ u which are 
singular in the limit g M — > 0; we will come back to this point in the next section. We thus 
can make for class B the following ansatz for the q dependence: 

T%(P, q, q') = a*" + b^q p + d» a > v q p q + ■■■. (3.6) 

The coefficients a ,i,v , • • ■ are 4x4 matrices which have to be constructed from the 16 Dirac 
matrices, the metric tensor g a/3 , the completely antisymmetric Levi-Civita pseudo tensor 
e a/37<s , and the remaining independent variables P a , and q' a . The form of these coefficients 
will be constrained by Lorentz invariance, gauge invariance, crossing symmetry and the 
discrete symmetries, C, P, and T. Contracting this ansatz with g M , the condition for the 
class B operator, Eq. (|3.5|), becomes 

q^ = a^q p + b^q p q p + • • • = 0. (3.7) 

Multiple partial differentiation of Eq. ( |3.7|) with respect to q a results in the following condi- 
tions for the coefficients, 

= 0, few," + = 0; c wa ' v = 0, • • • • (3.8) 

6 perm. 

Since current conservation is expected to hold for arbitrary q, the technique of partial differ- 
entiation to obtain conditions for the coefficients can easily be extended to obtain constraints 
for higher-order coefficients. The simple constraints of Eq. (|3.8|) are based on the fact that 
P a , q a , and q' a are independent variables; an implicit dependence would make matters more 
complicated. 

From a^ ,u = we can already conclude that the operator of class B contains no contri- 
butions which involve powers of q' only, without powers of q. Taking Eq. ( |3.8j ) into account, 
we now expand with respect to q', 

r%{P, q, q') = B^q p + B^ a q p q' a + ■■■ + C^q p q a + C^ a q p q a q' a + ■■■. (3.9) 

If we apply crossing symmetry to Eq. (|3.9| ) we find, as expected, that indeed all terms vanish 
that involve powers of q' only. Thus we find for the leading term of the model-dependent 
class B, 

T%(P, q, q') = B^ ua {P)q p q' a + 0{qqq', qq'q'), B m ™ = -B^ va = B ua ^ p , (3.10) 



In the following we omit the subscript "e//" . 
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where the conditions for B^ p,ua result from gauge invariance and crossing symmetry, re- 
spectively. To be specific, after imposing the constraints of Eq. (|3.10|) and of the discrete 
symmetries one obtains three possible structures for B^ p ' ua on the operator level: 

+i(g pu P p P a + g pa P»P v - gPvpup* - g pa pPp»)f 2 (P 2 ) 
+^ a lMP 2 ), 60123 = 1- (3.H) 

In summary, we have shown that on the operator level 

• the terms of order 0{q~ l ), 0{q'~ l ), 0(1), 0(q a ) and 0(q' a ) are contained in T^eft, Eq. 
(p.4j). They are therefore determined model independently through on-shell quantities; 



• model-dependent terms first appear in the order O^qaq'a); 

• all operators which contain either only powers of q or only powers of q' can entirely be 
obtained from this T^ e //- 

In the next section we will discuss the implications of these findings for the on-shell VCS 
matrix element. 



B. Application 

We now want to apply the above result to the observable case where the nucleons are on 
mass shell, i.e., we consider the matrix element. For that purpose we first define the matrix 
element of T pv between positive-energy spinors as 

V^ f (P,q,q , )=u(p f ,s f )T^(P,q,q') U (p i ,s i ). (3.12) 

At this point we have to keep in mind that for the variables we use, P, q, and q', the on-shell 
condition p 2 = pj = M 2 is equivalent to P ■ (q — q') =0, and P 2 + (q — q') 2 = AM 2 . In other 
words, the four-momenta chosen for the description of the off-shell Green's function will no 
longer be independent for the on-shell invariant amplitude. In particular, the distinction 
between powers of q only or, respectively, of q' only in T pu is not valid anymore for the 
matrix element since q ■ P = q' ■ P. 

Let us consider as an application the case where the initial photon is virtual and spacelike 
and the final photon is real, j*(q, e)+p(Pi, s») — ► j(q', e ')+p(Pf, Sf). The following discussion 
does not include the Bethe-Heitler terms of the physical process p(e, e'p)j, where the real 
photon is radiated by the initial or final electron, since these terms are not part of T pv . For 
the final photon the Lorentz condition q' ■ e' = is automatically satisfied, and we are free to 
choose, in addition, the Coulomb gauge e ,fl = (0, e') which implies q ' ■ e ' — 0. We write the 
invariant amplitude in the convention of |28[] as M. = — ie 2 e fJi M fl , where e M = eu^ u ujq 2 is 



the polarization vector of the virtual photon. With a suitable choice for the reference frame, 
<f — (<7o> 0, 0, \q\), the Lorentz condition q ■ e = and current conservation, g M M M = 0, may 
be used to reexpress the invariant amplitude as |p9fl 

M = ie 2 • M T + ^e 2 M, V (3.13) 
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Note that as go - > 0, both e z and M z tend to zero such that M. in Eq. Q3.13| ) remains finite. 
Making use that we now know the general structure of the first undetermined operator, Eq. 
( p. 11| ), we can now consider the matrix element up to second order in q or q'. This will 
enable us to explicitly show how far the amplitude is determined and where the first model- 
dependent terms enter. After a reduction from Dirac spinors to Pauli spinors the transverse 
and longitudinal parts of Eq. (|3.13|) may be expressed in terms of 8 and 4 independent 
structures, respectively 

e T ■ M T = e'* ■ e T A 1 + ia ■ (e'* x e T )A 2 + (<?' x e'*) ■ (q x e T )A 3 + ia ■ (q' x e'*) x(?x e T )A A 
+iq ■ e'*a ■ (q x e T )A 5 + iq ■ e T o ■ (q x e'*)A e + iq ■ e'*a ■ (q x ex)A 7 
+iq' ■ e T a ■ (q x e'*)A 8 , (3.14) 
e z M z = e z e'* ■ qA 9 + ie z e'* ■ qa ■ {q' x q)A w + ie z a ■ (q x e'*)A u + ie z a ■ {q' x e'*)A 12 . 

(3.15) 

The results for the functions A i in the CM system expanded up to 0(2), i.e. \q'\ 2 , \q'\\q\ 
and \q\ 2 , are shown in Tables | and [TT]. The derivation of the corresponding expressions is 
outlined in Appendix A. 

In the amplitudes A§ and Ag, terms of order l/|(f'| appear. We have not further expanded 
them in \q\ and kept, e.g., the on-shell form factors Gm and Ge- This was done mainly to 
stress that according to Low's theorem |J these divergent terms are already entirely fixed. 
They are due to soft radiation off external lines, with the intermediate line approaching the 
mass shell in the pole terms. They can be obtained from, e.g., the Born terms that contain 
the same on-shell information (see however the caveats in Sec. IV). In order to uniquely 
identify these singular contributions we have expanded all relevant expressions in terms of 
\q\ and \q'\. This is the reason why the argument of the form factors is Q 2 = g 2 ||gV|=o = 
—2M(Ei - M), where E { = V / AF + ?. 

Not only the irregular contribution, but all terms in the amplitudes up to terms linear 
in the photon momenta are uniquely determined through well-known properties of the free 
nucleon: its mass, charge and magnetic moment. These are the terms that make up the 
LET for VCS. Up to 0(2), the coefficients in addition also involve the electric mean square 
radius, which we know from electron-proton scattering, as well as the electric and magnetic 
polarizabilities which also enter in real Compton scattering. The latter are the coefficients 
of the first model-dependent terms in the expansion of the scattering matrix element. Their 
specific form, and thus the definitions of the polarizabilities, depend on the particular repre- 
sentation one chooses for the effective Lagrangian; the 0(2) results in the Tables are specific 
for the "canonical form" , which happens to be the standard choice for real Compton scatter- 
ing. Note that to the order considered A? = —Ag = Aiq. The amplitudes A\ and Ag contain 
the electric polarizability a of the proton whereas A3 involves the magnetic polarizability (3 
Q. In terms of the functions fi(P 2 ) of Eq. (|3.11| ) they are defined as 

a = e 2 (/i(4M 2 ) + 4M 2 / 2 (4M 2 )), = -e 2 /i(4M 2 ). (3.16) 

The function /s(4M 2 ) only contributes at order 0(3), since 75 connects upper and lower 
component of positive energy spinors which effectively leads to an additional power of |g| or 
\q'\ in the matrix element. 
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In Table |TXT| we also show the results for the transverse amplitudes in the CM frame for 
real Compton scattering. Since A 5 = —A 6 and A 7 = —A 8 the result effectively involves 6 
independent amplitudes as required by time reversal invariance 0. When comparing with 
other expressions in the literature flp|fl,|30| , one has to keep in mind that they usually are 
given in the laboratory frame. This observation accounts, for example, for the difference of 
the LET for A% and A% of real Compton scattering in the lab and the CM frame. Note also 
that the l/\q'\ singularity disappears in the real Compton scattering limit, since ui = \q\ = 
\q'\ in this case. 

In conclusion, the low-energy behavior of the VCS matrix element for e~ +p — > e~+p+7, 
expanded up to order 0(2) in \q\ and \q'\, contains, in addition to the structure coefficients 
that enter into real Compton scattering, also the electric mean square radius and the electric 
and magnetic Sachs form factors in the spacelike region which all can be obtained from 



electron scattering off the proton. For the reaction 7 + p — > p + e + e~ |13[ , the analogous 
information for timelike momentum transfers is needed. However, here one has to keep in 
mind that this information is not directly accessible for < q 2 < 4M 2 . 



IV. THE BORN TERMS 

From the above discussion in a particular representation one might be tempted to con- 
clude that the leading terms for VCS - on the operator level or for the amplitude - are in 
general simply given by "Born terms". However, some caveats are in order. First, one has 
to keep in mind that the low-energy behavior was obtained from considerations involving 
the most general ansatz for the truncated four-point Green's function. All terms one can 
think of are included into class A and class B. In, e.g., the derivation of || of the LET for 
Compton scattering, where no special representation is chosen, it is clearly shown that both 
class A and class B terms are needed. Implicit in all derivations is of course the assumption 
that a description in terms of observable asymptotic hadronic degrees of freedom is sufficient 
and complete at low energies. Even though subnucleonic degrees of freedom, quarks and 
gluons, are ultimately the origin of the structure of the nucleon, it was shown [0,0 that 
an effective field-theory approach P^-|57| in terms of hadrons is meaningful at low energies, 



thus allowing a classification into class A and class B. 

Secondly, there is an ambiguity concerning what exactly is meant by "Born terms" , once 
phenomenological form factors are introduced and the result is not obtained from a micro- 
scopic Lagrangian. We will illustrate this ambiguity by considering different representations 
of the photon-nucleon vertex. All representations contain the same information concern- 
ing the electromagnetic structure of the on-shell nucleon, as obtained in electron-nucleon 
scattering, but differ in the half-off-shell situation encountered in the s- and u-channel pole 
terms of Compton scattering. This difference is, of course, accompanied by different class B 
terms such that the total result is the same. 

To explain the above points in more detail, we will first reconsider the most general 
expression for class A and, without going to a special representation, identify those terms 
which contain the irregular contribution for g M — > or q'^ —>■ 0. We find that these contribu- 
tions can be expressed in terms of on-shell quantities, in our case the Dirac and Pauli form 
factors F\ and F 2 . We then show that the use of on-shell equivalent electromagnetic vertices 
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gives rise to the same results for the VCS matrix element as far as the irregular terms are 
concerned, but not every choice will result in "Born terms" which are gauge invariant. 



A. Irregular Contribution to the VCS matrix element 

When calculating VJ^ , the irregular contribution originates from the singularities of the 
propagators S(pi + q) and S(pi — q') in the s and u channel, respectively. To be specific, 
below pion-production threshold the renormalized, full propagator can be written as 

S(p + q) = Sp(p + q) + regular terms, (p + q) 2 < (M + m n ) 2 , (4.1) 

where Sf{p) denotes the free propagator of a nucleon with mass M. By "regular terms" 
we mean terms which have a well-defined, non-singular limit for q^ — > 0. Thus, as long as 
we are only interested in the irregular terms, we can simply replace the full, renormalized 
propagator by the free Feynman propagator, Sp. 

The most general form of the irreducible, electromagnetic vertex of the nucleon can be 



expressed in terms of 12 operators and associated form functions p8| , p9| . These functions 
depend on three scalar variables, e.g., the squared momentum transfer and the invariant 
masses of the initial and final nucleon lines. A convenient parametrization of r M (p/,pj) is 
given by 



r*(p/,ft)= E UPf)UFf + l ^Ff + ^Ff)Mp t ),F^ = F^(q 2 y f ,p 2 



a, /3= 



i ) i 



(4.2) 



where q = Pf — Pi, and A±(p) = (M ±tf)/2M. We have chosen a form which differs 
slightly from the convention of |39j in the definition of the projection operators and the 
normalization of the F 3 form functions. We only need the following on-shell properties of 
the form functions: 

F++(g 2 ,M 2 ,M 2 ) = F 1 (q 2 ), F++ (q 2 , M 2 , M 2 ) = F 2 (q 2 ), F 3 ++ (q 2 , M 2 , M 2 ) = F 3 (q 2 ) = 0, 

(4.3) 

where Fx and F 2 are the standard Dirac and Pauli form factors and F 3 vanishes because of 
time-reversal invariance. One can also show that F 3 (q 2 ) vanishes due to current conservation. 
We now systematically isolate the irregular part of, e.g., the s-channel pole diagram,^ 

v Zs = u(Pf) TU (Pf>Pf + q') iS (Pi + q) Tfl (Pi + q,Pi)u(Pi) 

w u(pf)T u (p f ,p f + q')iS F (pi + q)T^(pi + q,Pi)u(pi), 

where we made use of Eq. ( f4.1| ), and where the symbol ~ denotes equality up to regular terms. 
We now insert Eq. ( |4.2|) for T M (pj + q,Pi) and use A_(pi)u(pi) = and A + (pi)u(pi) = u{pi) 
to obtain 



2 In the following we omit the indices Sj and Sf. 
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VZ « u(p f )T v ip f ,p f + q')iS F { Pi + q) ( A+fe + q){YFf + {q\ s, M 2 ) + ■ • •) 

+ A_fe + g)( 7 M Ff + (g 2 , s, M 2 ) + • • .))u(p»), 
where s = (p* + q) 2 . Since Sjr(Pt + Q)-^—iPi + ?) results in a regular term, we have 

We now expand the form functions around so = M 2 and note that in the higher-order terms 
the powers of (s — M 2 ) cancel the denominator, s — M 2 , of the Feynman propagator and 
thus give rise to regular terms, 

VZ ~ u(p f )r»(p f ,p f + q')iS F { Pi + g)A + fe + qX-fFitf) + i-fiff 



where we made use of Eq. (|4.3|). Using A + ( Pi + q) — 1 — A_(pj + q) and then repeating the 
same procedure for T u (pj,pf + q') we finally obtain 

VZ ~ «(P/, Sf)r Fl>F2 (-q')iS F ( Pi + q)T FuF2 (q)u( Pl , Si ), (4-4) 
where we introduced the abbreviation 

rSi^(ff) = 7^1 (g 2 ) + i^f^iv 2 )- (4-5) 

The procedure for the u-channel part, V^, is completely analogous and we obtain for the 
sum of s- and u-channel contributions 

vz » n(p / )(r^ 1)F2 (- g / )i^(p« + ^^(g) + r^C?).^ - q ')r FuF2 (-gO)«(w) 

= (4.6) 

which only involves on-shell quantities, the Dirac and Pauli form factors and the nucleon 
mass. Explicit calculation, including the use of the Dirac equation, shows that VZ of Eq. 
(4.6) is, in fact, identical with evaluating Eq. (3.4) between on-shell spinors. With the 
help of either Eq. Q2.12p or by straightforward calculation it can easily be shown that Eq. 
( |4.6|) is gauge invariant. This is a special feature when working with these particular vertex 
operators and was essential for the derivation in |TJ|]. It is quite unexpected, since the 
electromagnetic vertex, Eq. (|4.5|) , and the nucleon propagator in Eq. (|4.6|) do not satisfy the 
Ward-Takahashi identity (except at the real photon point, q 2 = 0): 

q,r FuF2 ( Pf , Pt ) = Uj-tfFtf) ± s F \ P} ) - s F \ Pi ). (4.7) 



B. Different representations of the on-shell vertex 

We now turn to "Born-term" calculations involving other representations of the nucleon 
current operator, i.e., other ways to introduce the model-independent information about the 
electromagnetic structure of the free nucleon into the calculation of the lowest-order terms. 



Two commonly used alternative ways to parametrize the nucleon current operator are [40 
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r& lA (p/,ft) = 7^i(? 2 ) - ^^(g 2 ), (4.9) 

where P = Pi + pf and 

G E = F 1 + ^F 2 , G m = H 1 =F 1 + F 2 , H 2 = F 2 . (4.10) 

Given Eqs. ( |4.10| ), it is straightforward to show the equivalence for the free proton current, 
the matrix elements of Eqs. (|3.3|) , ( |4.5|) , (|4.8| ) and ( f4.9|) between free positive-energy spinors. 
On the other hand, the current operators in Eqs. ( ^4.5| ), ( fL8| ) and ( |4.9| ) do not satisfy the 
Ward-Takahashi identity when used in conjunction with free propagators, not even at the 
real-photon point (recall Eq. ( [4.7] )): 

%^G E ,G M (Pf,Pi) = (l " 4^2) 1 ^T G ^ * S ^ - ( 411 ) 

2 2 

<l^ Hl , H MPi) = (^f-^)Hi(q 2 ) - P -^T H ^ 2 ) * - S- F \ Pl ). (4.12) 

The "Born terms" calculated with these electromagnetic vertices and free nucleon prop- 
agators are 

v x = u{Pf,s f )(T v x (p } ,p f + q')iS F (pi + q)T^(pi + q,pi) 

+ T x(Pf>Pf - q)iS F (pi - q')T u x (pi - q',pi))u(pi, s<), (4.13) 



where X denotes either the vertex of Eq. (|4.8|) in terms of Ge, Gm or the vertex of Eq. ( |4.9| ) 
involving Hi,H 2 , respectively. These "Born terms" are by construction crossing symmetric 
but in both cases not gauge invariant, 



% V G E ,G U =i(l- j|p J G E (q 2 )u(p f , s f )(T u GEtGM (p f , Pi + q)- r% BtGu (pf ~ q,Pi) 

+r GB ,G M (Pf,Pi + q)^j + ^f GE ,G M {pf - q,Pi)H Pi , «,), (4.w) 

q» V HuH 2 = ^H^q 2 ) - H 2 (q 2 ))u(p f , s f )(T u HuH2 (p f , Pi + q) - T v HuH2 {p f - q,Pi))u(pi, S;) 

- zH 2 (q 2 )u(p f , s f )(r» HitH2 (p f , Pi + q)^L + JL r ^ ^{p } - q,p f ))u( Pi , s t ). (4.15) 



With the help of Eqs. (|4.8|) and ( |4.9|) it is easy to see that in each case the non-vanishing 
divergence is of order q. The fact that Eqs. (|4.14|) and ( 4.15|) do not contain terms of the 
type a ■ q/b ■ q suggests that in general the different "Born terms", based on equivalent 
parametrizations of the on-shell currents, will differ by regular terms, only. It is of course 
conceivable that the amplitude could differ by irregular terms, which are separately gauge 
invariant. Thus the above argument is not stringent and can only serve as a motivation for 
the following claim which is essentially equivalent to Low's theorem applied to the particular 
case of VCS: Any "Born-term" calculation involving electromagnetic current operators which 
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correctly reproduce the on-shell electromagnetic current of the nucleon will yield the same 
irregular contribution to the VCS matrix element. The key to the proof of this statement is 
the fact that any current operator which transforms as a Lorentz four- vector can be brought 
into the form of Eq. ( |4.2|) . On-shell equivalence then amounts to the constraint that all 
operators have the same on-shell limit of the Ff + form functions. In general, no statement 
can be made for either the other form functions or off-shell kinematics.0 However, as we have 
seen above, the irregular contribution of class A, and thus of the total VCS matrix element, 
only involves the on-shell information contained in F 1 ++ (g 2 , M 2 , M 2 ) and F 2 ++ (q 2 , M 2 , M 2 ). 
Any information beyond this will give rise to regular terms. Thus the above statement is 
true. 

To illustrate the above, we bring for example the current operators of Eq. (|4~5| ) and ( |4.9|) . 
involving F\ and F2 or Hi and H 2 , respectively, into the general form of Eq. ( [4.2| ) . By using 
1 = A + (p) + A_(p), the result for the commonly used form of Eq. (|4.5|) is given by 

Ff(q 2 ,p}, P t)=F 1 ( q 2 ), Ff( q 2 ,p 2 ,p 2 ) = F 2 (q 2 ), Ff(q 2 jP 2 ,p 2 )=0, a,(3 = +,-. 

(4.16) 

For the vertex given in Eq. ( |4.9| ), we use {7^,7^} = 2g fll/ and momentum conservation at 
the vertex to rewrite 

2M 2M 1 2M ' 1 ' ' 

By inserting appropriate projection operators in the form p 1 = M(A + (p) — A_(p)) and as 
above, the vertex of Eq. (|49|) can be expressed as 

^ Hl>H2 (Pf,Pi) = A +(P/) (rma- 2 ) - H 2 {q 2 )) + f-^H 2 { q 2 )) A + (pA 
+ A_(p,) ( 7 ^i(<7 2 ) + i W^(g 2 )) A + fe) 

2^ , .^q v 



A+(P/) ( YH.iq 2 ) + *^H 2 (q 2 ) ) A_fe) 



+ A_(p f ) ^(H.iq 2 ) + H 2 (q 2 )) + *^fH 2 (q 2 ) j A_(ft). (4.18) 

Of course, Eq. ( [4. 18|) contains the same on-shell information as Eq. ( |4.5|) , since the form 
factors satisfy Eq. (|4.10|) . On the other hand, the expressions for all the other F^ 13 form 
functions differ. This is why the two "Born term" calculations based on these two vertices 
differ with respect to regular terms. It is straightforward to extend the same considerations 
to the vertex involving Ge and Gm- 

In conclusion, we have shown that a calculation based on only the "Born terms", built 
from any of the many possible on-shell equivalent vertices and free nucleon propagators, 
yields the same results for the irregular terms as the LET. Thus these "Born terms" will 



3 Further conditions on the form functions can be derived from the Ward-Takahashi identity and 



discrete symmetry requirements [ 38,39 1. 



16 



differ among each other through regular terms. Furthermore, such "Born terms" are, in 
general, not gauge invariant; an exception is the commonly used form involving the Dirac 
and Pauli form factors F\ and F 2 . "Generalized Born terms" which are made gauge invariant 
by hand through an ad hoc prescription also differ by regular terms. 

Important starting point for the derivation of the LET are the irregular terms. It is 
thus also possible to split the total VCS amplitude into "Born terms" plus "rest" , instead of 
class A and B amplitudes, to arrive at the same result for the LET, i.e. up to and including 
terms linear in the photon three-momenta. In general, this result will have contributions 
from "Born terms" and the "rest" amplitude. If one uses a "generalized Born amplitude", 
all the terms appearing in the LET are due to the expansion of the Born amplitude. It is a 
well-known feature of soft-photon theorems that they cannot make statements about terms 
which are separately gauge invariant [ jl- j^. One has to keep this in mind when discussing 
the structure-dependent higher-order terms of VCS, i.e., one needs to specify which Born 
or class- A terms have been separated. For example, in |TJJ the "Born terms" involving 



Fi and F 2 where separated since they provide without any further manipulation a gauge- 
invariant amplitude. Then the residual part with respect to these particular "Born terms" 
was parametrized in terms of generalized polarizabilities. A natural question to ask is what 
would have happened had one separated a different choice of "generalized Born terms" and 
defined generalized polarizabilities in an analogous fashion with respect to the corresponding 
residual amplitude. Obviously one would, in general, have found different numerical values 
for the new generalized polarizabilities in order to obtain the same total result. 



V. CONCLUSIONS 

In studying the structure of composite strongly interacting systems the electromagnetic 
interaction has been the traditional and precise tool of investigation. In scattering of elec- 
trons from a nucleon, our knowledge is restricted to two form factors that we can extract 
from experiments. Even though we have not yet been able to fully explain this information 
on the basis of QCD, it is important to look for other observables allowing us to test ap- 
proximations to the exact QCD solution and effective, QCD inspired models. Such effective 
models are expected to work especially at low energies. The electron accelerators now make 
it possible to study virtual Compton scattering, which is clearly more powerful in probing 
the nucleon than the scattering of real photons. In analyzing Compton scattering it is im- 
portant to know how much of the prediction is not a true test of a model, but fixed due 
to general principles. These model-independent predictions for virtual Compton scattering 
were the main topic of our discussion. 

The interest in virtual Compton scattering has also been due to another aspect: When 
studying reactions on a nucleus, such as (e,e'p), the nucleon interacting with the electro- 
magnetic probe is necessarily off its mass shell. We have no model-independent information 
for the behavior of such a nucleon and any conclusion about genuine medium modifications 
must be based on firm theoretical ground how to deal with a single nucleon under these kine- 
matical circumstances. In fact, such a discussion depends very much on what one chooses 
as an interpolating field for the intermediate, not observed nucleon. This clearly makes the 
"off-shell behavior" of the nucleon representation dependent and unobservable. We discussed 
how certain features of the off-shell electromagnetic vertex of the nucleon can be shifted into 
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irreducible, reaction-specific terms for the reaction amplitude. Two-step reactions on a free 
nucleon - like (virtual) Compton scattering - allow us to test many aspects of dealing with 
an intermediate, off-shell nucleon under simpler circumstances, without complications from, 
e.g., exchange currents or final state interactions. Understanding these aspects on the single 
nucleon level would seem a prerequisite before any exotic claims can be made for nuclear 
reactions. 

We have studied the virtual Compton scattering first on the operator level. Using the 
requirement of gauge invariance, as expressed by the Ward-Takahashi identity, we derived 
constraints for the operator that determine terms up to and including linear in the four- 
momenta q and q' . Also, we showed that on the operator level terms involving terms de- 
pending only on q or only on q' are determined model independently in terms of on-shell 
properties of the nucleon. 

To obtain these results, we used the method of Gell-Mann and Goldberger, by splitting 
the contributions into general pole terms (class A) and the one-particle irreducible two- 
photon contributions (class B). We calculated class A below pion-production threshold in 
the framework of a specific representation for the most general effective Lagrangian compat- 
ible with Lorentz invariance, gauge invariance and discrete symmetries. This approach was 
introduced in M as a method for writing the general structure of the Compton scattering 
amplitude in a way that allows best to discuss its low-energy behavior. In this connection, 
we also showed the origin for a commonly used form of the electromagetic vertex of the 
nucleon and stated the consistency conditions for its use. 

After discussing the leading terms of the VCS operator, we considered the matrix element 
for 7*p — > 7p in the photon-nucleon CM frame. We found that the VCS amplitude up to 
and including terms linear in the initial and final photon three-momentum can be expressed 
in terms of information one can obtain from electron-proton scattering. This is the result 
analogous to the LET for the real Compton scattering amplitude. As we also showed, the 
next order - terms involving |<f' | 2 , |<f' and \q\ 2 - is also completely specified but now 
requires in addition also the electromagnetic polarizabilities a and (3 encountered in real 
Compton scattering. In other words, new structure-dependent information can only appear 
at order three or higher in the three-momenta. Our results concern the expansion in terms of 
powers of both the initial and final photon momentum. This allowed us to determine more 
terms than in ]TJJ , where only the leading terms in the final momentum were concerned. On 
the other hand, by expanding in both momenta, the range of applicability is smaller since 
both kinematical variables should be small. 

We then considered different commonly used methods to include the on-shell information 
contained in the electromagnetic form factors in a "Born-term" calculation of the VCS matrix 
element. The fact that the "Born terms" calculated with F\ and F2 are gauge invariant is not 
trivial, since the vertex and free propagator do not satisfy the Ward-Takahashi identity. We 
explained why different on-shell equivalent forms for the electromagnetic vertex operator 
lead to the same irregular contribution in the VCS matrix element. We emphasized the 
importance of stating with respect to which pole terms the structure-dependent terms are 
defined. 

Using only gauge invariance, Lorentz invariance, crossing symmetry, and the discrete 
symmetries, we were able to make statements about the low-energy behavior up to 0(2). 
Further conclusions can be reached by also taking into account the constraints imposed by 
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chiral symmetry. This would most naturally be done in the framework of Chiral Perturbation 
Theory. In particular, predictions for the higher-order terms could be obtained. 
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APPENDIX A: 

In this appendix we outline the calculation of the transverse and longitudinal functions 
Ai of Eq. ( |3.14| ) and ( |3.15| ), respectively. For that purpose we split into its contributions 
from classes A and B, = + M^. If we introduce 

^(-90 = ^(1 + 2^^0 (Al) 



for the vertex involving the final real photon, the contribution of class A (see Eq. Q3.4Q ) 
reads: 



M% = u(p f ,8 f ) (F(-q')S F ( Pi + q)T' ff (q) + T» ff (q)S F ( Pi - q')F(-q')) u( Pi , Si ), (A2) 
where T^j(q) is defined in Eq. (^.3|). Applying the Dirac equation, can be written as 

= U{JP ^ ] 2 (7=Mi + ^W 2 V ^ q) + {1 + K) s-M 2 + u-M 2 



M s-M 2 u-M 2 + 2M ( t T *ff^ + r effW ) 



(A3) 



where s = (pj + q') 2 and u = (pi — q') 2 . Similarly, using Eq. (|3.10|) and ( |3.11|) , can be 
written as 

= u(p f , 8f ) ((e'^q ■ q' - g"Y* • q)h(P 2 ) 

+ (e'* M P -qP-q + P^e* ■ Pq ■ q' - PV* ■ qP ■ q' - g'V* ■ PP ■ q)f 2 {P 2 ) 
- ie ^e':q P q a ^h{P 2 ) + 0(3)) u(p t , s t ), (A4) 



where 0(3) denotes terms of order 3 in q or q'. 
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The following considerations will be carried out in the center-of-mass (CM) frame, where 
the four-momenta are given by g M = (g , 0,0, |g|), p% = (E g , —q), q ,fl = (\q'\,q'), and = 
(E q >, —q'). From energy conservation, q + E q = \q'\ + E g i, we infer that we may choose 
\q\, \q'\, and z = q ■ <f as a set of independent variables. In terms of the CM variables the 
denominators of Eq. (|A3|) are proportional to \q'\, 



M 2 = 2\q'\{E q , + |g'|), u — M 2 = -2\q'\(E q + z\q\), (A5) 



and thus can be written as 



M A = ^f + b^q,q'). (A6) 



The functions a^{q, q') and b^{q, q 1 ) are regular with respect to \q'\ and |gj, and are given by 
aP(q,<t) = K(q)u(-q',s f )T^ f (q)u(-q, Sl ), (A7) 
fc"(g, q') = u(-q', a,) ({1 + k) 



2{E q , + \q>\) 2(E q + z\q\) 



-^W 2§t|) + 2^'* r ^ + ) «>' (A8) 



where we introduced = (1, <f), and 

K(q) = -Tr^r- (A9) 
£ g + z|g| 

In order to obtain Eq. ( |A9| ) we explicitly made use of the Coulomb gauge for the final photon, 
namely, when using pf ■ e'* = and Pi ■ e'* = q ■ e'*. 

We will now derive from Eqs. (|A6|) - (|A9|) the expansion of in powers of |g| and \q'\. 
According to Eq. Q3.13| ) it is sufficient to treat the space components of M^. Since the struc- 
ture coefficients appear at 0(qq') and higher in the operator, we expand the contributions 
up to and including \q'\ 2 , \q'\\q\ and |g| 2 . When expanding the electromagnetic vertex we 
make use of the following relations 

<z" = Q" + |gV° + ---, 

4 = <#+\q'\io + -~, 

q 2 = Q 2 + 2Q \q'\ + ---, (A10) 

where = g M |,f=o = {M — Ei,q). Furthermore, for a we need the expansion of the form 
factors around Q 2 : 

F^q 2 ) = FiM) + 2Q \q'\Fi 2 (Q 2 ) + • • • , (All) 



with the notation F'(x) = dF/dx. Using the definition of T^^(q) of Eq. fl3.3|), it is straight- 
forward to obtain 

f e //(g) = f 00 + f 10 |g'| + f 01 |g1 + • • ■ , (A12) 
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where 



Too = 7> r 10 



AM 

470^(0), f 02 



[7, 7o] 



01 



4M 



[7,7-9], 



20 



^'(0)7 



8M 2 



[7, 7o] 



K 



(-7 + 7-^)^(0) + ^[7,7o] 



(A13) 



Since the dependence on the momenta q and q' is also contained in the initial and final 
nucleon spinors, respectively, we expand them to the required order: 



u(-q) 



:i + + iSr + • • -H0), u(-q') = tZ(0)(l + ■ 



1 9 



/|2 



2M 8M 2 



2M 8M 2 



+ 



(A14) 



where from now on we suppress the spin indices. Finally, the expansion of the energy 
denominators reads 



\Q 



f\2 



Eg, + \q'\ M M 2 2M 3 



+ 



1 



z\q\ , (2^ 2 -l)|gt 



E q + z\q\ M M 2 



2M 3 



+ •••. (A15) 



Let us first consider a(q, q') which we expand according to 



a(q,q') = a 00 (Q) + |<fVio + \q'\ a 20 + |<f'||<7|aii + \q'\ a 30 + \q'\ \q\a 2 \ + \q'\\q\ a 12 + 



Using the relations of Eqs. (|A10| ) - ( |A15| ) we obtain: 
a 00 (Q) = K(q)u(0)f eff (Q)u(-q), 



021 



Ol2 



q ■ e 
M~ 

til 
M 



-u(0) 



-u(0) 



r o + r 10 U(o), 



7-9' 
v 2M 

^r 00 + l^r 10 + r 20 



u(0), 



(A16) 

(A17) 
(A18) 

(A19) 



M 



M K 2M 



2M 



+fn + (^foo + f 10 ) 7 ' f/ 



2M 



2M 



u(0), 



I>j0 



0, J = 1,2, 



(A20) 
(A21) 



The last equation follows from K(0) = 0. The function Soo(Q) is completely determined in 
terms of the electromagnetic form factors Fi^iQ 2 ) (or Ge,m(Q 2 ))- Note that in Eq. ( |A17| ) 
we keep all powers in \q\, since it will be multiplied with the l/|g'| singularity and there are 
no other terms which can generate such a singularity. 

In order to determine b(q,q'), we first expand b(q,q') in an analogous fashion to Eq. 



b(q,q') = b 00 + \q'\b w + |# i + |<f R20 + |?'||#ii + \q\ %2 + 



(A22) 



Using the building blocks of Eq. ( |A10| ) - ( |A15| ) we find for the coefficients b^: 
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boo = 2^«(0) ((1 + k)W ,f 00 ] + ti(O) = ^«(0)Xtt(0), (A23) 

&10 = 2^«(0) ((1 + «)([^V ,f w ] - ^'V f 00 ) + «{^,f M } + «(0) 

boi = ^u(O) ((1 + «)(^' , f Q1 ] + ±f oc/V ) - ^ • gfV + «{^ , foi} + X^fj u(0) 

&20 = 2^«(0) ((1 + «)(^' , fa,] - f 10 + ^-J*i' f oo) + f 20 } 

'^Y + -^x)u{0), (A26) 



2M 8M 2 " 

&n = 2^«(0) ((1 + K)(/y , f n ] - Foi + ^iVV ) - • gfV + f n } 

■ « f « - s p »^ + kW " f ° 2} + z M + 8^ x ) u(o) - (A28) 

The reduction of the above expression to Pauli space is straightforward but very tedious. 
The results are displayed in Tables [TV| and [Vj 

Finally, it is straightforward to obtain the expansion of M#: 

M B = fi(0) [\q?e'*^ + \?M(zf - q ■ e'*g')| + 0(3)\ «(0), (A29) 

where we have defined the (real) Compton polarizabilities as 

a = e 2 (A(4M 2 ) + 4M 2 f 2 {4M 2 )), j3 = -e 2 /i(4M 2 ). (A30) 

Due to the presence of the 75 matrix, the third function f 3 (4M 2 ) only contributes at 0(3) 
at the level of the matrix element. 
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TABLES 



TABLE I. 

Transverse functions of Eq. ( |3.14| ) in the CM frame. The functions are expanded in 
terms of \q'\ and |g| of the final real and initial virtual photon, respectively. iVj = d / E ^ M 
is the normalization factor of the initial spinor, where Ei = \J M 2 + \q\ 2 . Ge{q 2 ) = Fi(q 2 ) + 

^p-F 2 (g 2 ) and Ga/(<? 2 ) = Fi(q 2 ) + F 2 (q 2 ) are the electric and magnetic Sachs form factors, 
respectively. r% = 6G' E (0) = (0.74±0.02) fm 2 is the electric mean square radius JJT] and k = 
1.79 the anomalous magnetic moment of the proton. Q 11 is defined as |q-'|=o = (M — Ei, q), 
Q 2 = —IMiEi — M), and z = q' ■ q. at and (3 are the electric and magnetic Compton 
polarizabilities of the proton, respectively. 



At 


1 1 ( 1 I r % K a\ I-/I2 | ( 1 | r% z 2 (1+k)k\ 1 -, 2 
M ^ I 8M^ ^ 6M iM 3 ' F J ^ 1 ^ I 8AF ^ 6M AF t 4Af ;i J 1*1 
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2M^ 19 1 4M^I9 1 1 2A/^I9 1191 4M 3 191 


A 3 


-^l9l + (4^ + l)l^k1 + ^S^k1 2 
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-^l9l-^l9'll9l + ^kT 
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(Bi+zMH^+iW) \q'\ 1 4Af 3 191 


A 6 
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A 7 




A 8 


^fl9"ll9l 



TABLE II. 



Longitudinal functions A\ of Eq. fl3.15| ) in the CM frame. See caption of Table |. 



A 9 


N t G E (Q 2 ) \q\ 2 1 
(Ei+z\q\)(Ei+M) \q'\ M 




.u _ { _i + li_ « _ 

l y l I 8M^ ~ 6M HP 


S ) I9-I 2 + (jfa + A " l9T 


A 10 




An 




1+2k 
2M 2 


9 1 + 4M?|9 1 + AM :i 19 1 


9l + ^l9T 


A 12 


{1+k)z 
2M 2 


q'\- 


(1+k)kz, -./| 2 (1+k)(2^-1) 
4M 3 19 1 4M 3 





25 



TABLE III. Transverse functions A4 in the CM frame for both photons real: q 2 = q' 2 = 0, 
\q\ = \q'\ = to. 
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TABLE IV. 

Reduction of the coefficients Oy of Eqs. ( |A17| ) - ( |A20|) to Pauli space. For the definition of 
the corresponding Pauli spin operators see Eqs. ( |3.14| ) and ( 3.15| ). For further information 
see caption of Table |. 
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Reduction of the coefficients b 



TABLE V. 

of Eqs. ( |A23|) - (|A28[) to Pauli space. See captions of Tables 



and IV 
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